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Abstract 



A search for quark compositeness using dijet angular distributions from pp collisions 
at \fs = 7 TeV is presented. The search has been carried out using a data sample cor- 
responding to an integrated luminosity of 2.2 fb -1 , recorded by the CMS experiment 
at the LHC. Normalized dijet angular distributions have been measured for dijet in- 
variant masses from 0.4 TeV to above 3 TeV and compared with a variety of contact 
interaction models, including those which take into account the effects of next-to- 
leading-order QCD corrections. The data are found to be in agreement with the pre- 
dictions of perturbative QCD, and lower limits are obtained on the contact interaction 
scale, ranging from 7.5 up to 14.5 TeV at 95% confidence level. 



Submitted to the Journal of High Energy Physics 



*See Appendix A for the list of collaboration members 



1 



In theories of physics beyond the standard model, it has been proposed that quarks are com- 
posite particles and are bound states of more fundamental entities |HJ[2]|. Models of quark com- 
positeness may explain the number of quark generations, quark charges, and quark masses, 
which are not predicted in the standard model. A common signature of quark compositeness 
models is the appearance of new interactions between quark constituents at a characteristic 
scale A that is much larger than the quark masses. At energies well below A, these interactions 
can be approximated by a contact interaction (CI) characterized by a four-fermion coupling. In 
this Letter, flavor-diagonal color-singlet couplings between quarks are studied. These can be 
described by the effective Lagrangian HH3l 



L qq = XT [^(q L V'qL)(q L 7?/qL) + ?/RR(q R 7 F qR)(q R 7FqR) +2j/RL(q R 7 F qR)(q L T/<qL)] , 

where the subscripts L and R refer to the chiral projections of the quark fields and rjn, rj RR , and 
r/Ri can be 0, +1, or —1. The various combinations of rjn, ^rr, and rj R i correspond to different 
CI models. The following CI scenarios are investigated: 



A = A^ L for (i]LL,f]RR,T] RL ) = (±1,0,0), 
A = A| R for (r] LLr rj RR ,t] RL ) = (0,±1,0), 
A = A vv for (tin, rj RR , r] RL ) = (±1, ±1, ±1), 
A = A AA for {y]Li,y] RR ,ri RL ) = (±1, ±1,^=1), 
A = A^_ A) for (f]it,r] RR ,r] RL ) = (0,0, ±1). 

In pp collisions these models result in the same limits for A^ L and A RR , and at tree level for 

Ay V and A^ A as well as for A^ V _ A ^ and A7y_ A y 

High energy proton-proton collisions with large momentum transfers predominantly produce 
events containing two jets with high transverse momenta (dijets). Such events probe the scatter- 
ing partons at the shortest distance scales and provide a fundamental test of quantum chromo- 
dynamics (QCD). The angular distribution of these two jets with respect to the beam direction 
is directly sensitive to the underlying dynamics of the parton-parton scattering and does not 
strongly depend on the parton distribution functions (PDFs). Distributions of the polar scat- 
tering angle 6* in the parton-parton center-of-mass frame from QCD processes are peaked in 
the forward and backward directions, whereas contact interactions give rise to more isotropic 
distributions in 9* . 

Previous searches for quark compositeness at hadron colliders have been reported at the SppS 
by the UA1 JH collaboration, at the Tevatron by the DO [5 , 6 j and CDF collaborations, and 
at the Large Hadron Collider (LHC) by the ATLAS HJH and CMS |10HTT1 collaborations. The 
limits on quark compositeness at the LHC l[8l-[TT| have been reported only for a color- and 
isospin-singlet CI model, A^ L/RR , where A Z I X/RR (A^ L/RR ) corresponds to destructive (con- 
structive) interference between the CI and QCD terms. In this Letter, our previous searches are 
extended to higher CI scales using a data sample corresponding to an integrated luminosity of 
2.2 fb -1 at ^/s = 7TeV, exploring for the first time at the LHC a wide range of CI models. Also, 
this is the first use of a recent CI prediction that includes next-to-leading-order (NLO) QCD 
corrections 02). 
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In this analysis, the normalized dijet angular distributions, defined as (l/Cdijet) (dc di j et / d^dijet) 
where ^dijet = el yi_y2 l, are studied for several ranges of the dijet invariant mass Mjj. Here, 
1/1 and 1/2 are the rapidities of the two highest transverse momentum (px) jets, and they are 
related to the jet energy E and the projection of the jet momentum on the beam axis, p z , by 
y = g ln[(£ + p z )/ (E — p z )]. In the limit of massless scattering partons, ^dijet is related to 9* 
by Xdijet = (1 + | cos 9* | ) / (1 — | cos 9*\). The use of the variable ^dijet is motivated by the fact 
that dcr di j et / d^dijet is approximately uniform for QCD dijet processes, while CI models predict 
angular distributions that are strongly peaked at low values of ^dijet- 

The data for this analysis are collected with the Compact Muon Solenoid (CMS) detector at 
the LHC. The central feature of the CMS detector is a superconducting solenoid, 12.5 m long 
and with an internal diameter of 6 m, providing an axial field of 3.8 T. The field volume of the 
solenoid is instrumented with various particle detection systems. Charged particle trajectories 
are measured by a silicon pixel and strip tracker, covering < q> < 2n in azimuth and \rj\ < 2.5, 
where pseudorapidity rj = — ln[tan(0/2)] and 9 is the polar angle relative to the counterclock- 
wise proton beam direction with respect to the center of the detector. A lead-tungstate crystal 
electromagnetic calorimeter and a brass/ scintillator hadronic calorimeter surround the track- 
ing volume. A preshower detector made of silicon sensor planes and lead absorbers is installed 
in front of the electromagnetic calorimeter at 1.653 < \tj\ < 2.6. Outside the solenoid, muons 
are measured in gas-ionization detectors embedded in the steel return yoke. A more detailed 
description of the CMS detector can be found elsewhere |Tl3| . 

The CMS detector records events with a two-tiered trigger system consisting of a hardware- 
based Level-1 (LI) and a software-based High Level Trigger (HLT). In this study, single-jet 
triggers that reconstruct jets from calorimeter energy deposits at LI and HLT are used to select 
events based on different jet-px thresholds. Seven combinations of (LI, HLT) px thresholds (in 
GeV) are used to select events: (36, 60), (68, 80), (92, 110), (92, 150), (92, 190), (92, 240), and 
(92, 300). All except the highest-threshold jet trigger were prescaled during the 2011 run. The 
efficiency of each single-jet trigger is measured as a function of Mjj using events selected by a 
lower-threshold trigger. 

Jets are reconstructed offline using the anti-fcx clustering algorithm with a distance parameter 
R = 0.5 [14J. The four-vectors of particles reconstructed by the CMS particle-flow algorithm 
are used as input to the jet-clustering algorithm. The particle-flow algorithm IT5llT6]l combines 
information from all CMS subdetectors to provide a complete list of long-lived particles in 
the event. Reconstructed and identified particles include muons, electrons (with associated 
bremsstrahlung photons), photons (including conversions in the tracker volume), and charged 
and neutral hadrons. The jet energy scale is calibrated using measurements of energy balance in 
dijet and photon+jet events fF7\ . Extra energy clustered into jets from additional proton-proton 
interactions within the same bunch crossing (pileup) is taken into account on an event-by-event 
basis by a correction to the jet four-vectors. The average number of pileup interactions for the 
data sample used in this analysis has been estimated to be 5. 

Events with at least two reconstructed jets are selected, and the two highest-px jets are used to 
measure the dijet angular distributions for different ranges in Mjj. Events with spurious jets 
from noise and non-collision backgrounds are rejected by applying loose quality criteria to the 
jet properties [18] and requiring a reconstructed primary vertex within ±24 cm of the detector 
center along the beam line and within 2 cm of the detector center in the plane transverse to the 
beam ffl9|. The rapidities \y\\ and |i/2| of the two highest-px jets are restricted to be less than 
2.5 by selecting only events with Xdijet < 16 and |y b oost| < 1-11/ where y hoost = \{y x +yi). 
The lower limits of the Mjj ranges for the dijet angular distributions were chosen such that 
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the trigger efficiencies exceed 99%, and are given by the values 0.4,0.6,0.8, 1.0, 1.2, 1.5, 1.9,2.4, 
and 3.0 TeV. The data for the first five Mjj ranges are recorded using prescaled triggers and 
correspond to integrated luminosities of 0.77, 5.9, 32, 108, and 371 pb _1 , while the data for the 
mass ranges with Mjj > 1.5 TeV correspond to the full integrated luminosity of 2.2 fb -1 . The 
uncertainty on the integrated luminosities has been estimated to be 4.5% |20H2T1. The highest 
value of Mjj observed in this data sample is 3.9 TeV. 

The dijet angular distributions are corrected for migration effects due to the finite jet energy 
and position resolutions. The four-momenta, rapidities, and azimuthal angles of generated jets 
from Monte Carlo (MC) event simulations are varied within their measured resolutions ||T7|, 
and correction factors for each Mjj region are obtained from the ratio of the generated to the 
smeared ^dijet distributions. Unfolding correction factors are evaluated from two independent 
MC samples, pythia 6.422 (221 with tune D6T [231 and HERWIG++ 2.4.2 ||24| with tune 2.3, 
and the average of these corrections is applied to the data. The size of the correction factors 
varies from less than 1.3% in the lowest Mjj range to less than 10% in the highest Mjj range. 
The associated systematic uncertainties are taken as the maximum differences between the un- 
folding corrections obtained from four independent MC samples, HERWIG++ tune 2.3, PYTHIA6 
tunes D6T and Z2 (the Z2 tune is identical to the Zl tune [23J except that Z2 uses the CTEQ6L 
PDF |25j), and PYTHIA8 [26J tune 4C [27], and the nominal correction factors. These uncertain- 
ties range from less than 0.2% at low Mjj to less than 4.9% at high Mjj. A systematic uncer- 
tainty from using a parameterized model to simulate the finite jet pj and position resolutions 
to determine the unfolding correction factors is estimated by comparing the smeared ^dijet dis- 
tributions to the ones from a detailed simulation of the CMS detector using GEANT4 [28|. This 
uncertainty is found to be less than 1.3% (2.0%) in the lowest (highest) Mu range and is added 
in quadrature to the unfolding uncertainties. 

The dijet angular distributions are normalized to the integrated dijet cross sections in each Mjj 
range and are relatively insensitive to many systematic effects. For example, they show little 
dependence on the overall jet-energy scale and are independent of the luminosity uncertainty. 
However, they are sensitive to the rapidity dependence of the jet energy calibration and to the 
jet px resolution. For the phase space in pj and t] of the jets in this analysis, the jet energy 
scale uncertainties vary between 2% and 3% and have a dependence on pseudorapidity of less 
than 1% per unit of rj |[17|. The uncertainty on the jet pj resolution is less than 10% [17J. The 
resulting uncertainty on the Xdijet distributions due to the jet energy calibration uncertainties is 
found to be less than 1.0% at low Mu and less than 0.3% at high Mjj over all ^dijet bins, while 
the maximum uncertainty due to the jet pj resolution uncertainty varies from 0.2% at low Mjj 
to 0.6% at high Mu. In addition, uncertainties on the tails of the jet pj resolutions [17] result in 
systematic uncertainties on the Xdijet distributions ranging from less than 0.5% at low Mjj to less 
than 4.6% at high Mjj. The effect of pileup was investigated by dividing the data into low and 
high pileup samples based upon the vertex multiplicity, and comparing the ^dijet distributions 
from each sample. No significant effect was observed. The total systematic uncertainty on the 
/Cdijet distributions, calculated as the quadratic sum of the contributions due to the uncertainties 
in the jet energy calibration, the jet pj resolution, and the unfolding correction, is less than 1.7% 
for the lowest Mjj range and less than 7% for the highest Mjj range. A summary of the leading 
systematic uncertainties is provided in Table [I] 

Predictions at NLO in perturbative QCD are made for the dijet angular distributions with NLO- 
JET++ 2.0.1 |29l in the FASTNLO framework version 1.4 |13"0"1 . The factorization (jiA and renor- 
malization (}i r ) scales are set to {pj), the mean pj of the two jets, and the PDFs are taken from 
the CTEQ6.6 set il3TH . Correction factors are applied to the predictions to account for non- 
perturbative effects due to hadronization and multiple parton interactions. These correction 
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Table 1: Summary of the leading experimental and theoretical uncertainties on the Xdi\et distri- 
butions. The maximum uncertainties over all Xdi]et bins for the lowest and highest M« ranges 
are given. The dominant experimental contribution is the statistical uncertainty while the dom- 
inant contribution to the NLO perturbative QCD uncertainty is the factorization and renormal- 
ization scale uncertainty. 



Source of Uncertainty 


0.4 < Mjj < 0.6 TeV 
(%) 


Mjj > 3 TeV 
(%) 


Jet energy scale 


1.0 


0.3 


Jet energy resolution 


0.2 


0.6 


Jet energy resolution tails 


0.5 


4.6 


Unfolding, MC modeling 


0.2 


4.9 


Unfolding, detector simulation 


1.3 


2.0 


Total experimental systematic uncertainty 


1.7 


7.0 


Statistical uncertainty 


2.5 


31.6 


]A r and y. j scales 


5.6 


14.9 


PDF (CTEQ6.6) 


0.5 


0.7 


Non-perturbative corrections 


1.7 


1.1 


Total theoretical systematic uncertainty 


5.9 


15.0 



factors are used to correct the parton QCD calculations to the particle level, and they are de- 
termined by the average of the corrections estimated using PYTHIA6 tune D6T and HERWIG++ 
tune 2.3. This uncertainty is found to encompass alternative choices of MC tunes, PYTHIA6 
tune Z2 or PYTHIA8 tune 4C, and is estimated to be less than 1.7% (1.1%) at low (high) Mjj. 

The dominant source of uncertainty on the QCD predictions is due to the choices of the ji t 
and ji r scales. The uncertainty is evaluated following the proposal in Ref. [32] by varying 
the default choice of scales in the following 6 combinations: (}if/ {pj}, fi r / (Pt)) = (1/2, 1/2), 
(1/2,1), (1,1/2), (2,2), (2,1) and (1,2). These scale variations modify the predictions of the 
normalized ^dijet distributions by less than 5.6% (15%) at low (high) Mjj. The uncertainty due 
to the choice of PDFs is determined from the 22 uncertainty eigenvectors of CTEQ6.6 using the 
procedure described in Ref. [31 j, and is found to be less than 0.5% at low Mjj and less than 
0.7% at high Mjj. The leading systematic uncertainties on the theoretical predictions are listed 
in Tabled] 

The measured differential dijet angular distributions, corrected for instrumental effects and 
normalized to their respective integrals, are compared to QCD predictions in Fig.[I]for different 
Mjj ranges. Overall the theoretical predictions provide a good description of the data for all Mjj 
ranges. 

The measured dijet angular distributions are used to set limits on a variety of CI models. Only 
color-singlet models, which predict the largest deviations of the dijet angular distributions from 
the standard model, are considered. In fact, for the general case of a CI model containing both 
color-singlet and color-octet contributions, there are certain regions in the theory parameter 
space where the CI predictions for the dijet angular distributions become indistinguishable 
from the QCD predictions, because of interference between these contributions |Tl2| . 

In this analysis we present limits for a CI model that includes the exact NLO QCD corrections 
to dijet production induced by contact interactions |Tl2| , as well as limits extracted from various 
CI models implemented in PYTHIA8 [26]. In the latter case, the contributions of CI and QCD 
are calculated to leading order (LO). To take into account the NLO QCD corrections which 
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Figure 1: Normalized dijet angular distributions for | i/boost I < 1.11 in several M« ranges. For 
clarity, the distributions are shifted vertically by the additive amounts shown in parentheses in 
the figure. The vertical error bars represent the statistical and systematic uncertainties added 
in quadrature. The horizontal bars correspond to the ^dijet Drn width. The results are compared 
with the predictions of NLO QCD with CTEQ6.6 PDF (shaded band) and with predictions for 
QCD+CI from [12] at the CI scale A^ L/ , RR = 7 TeV (dashed histogram). Non-perturbative cor- 
rections due to hadronization and multiple parton interactions are applied to the predictions. 
The shaded band indicates the total uncertainty on the NLO QCD predictions due to ]i r and 
ji f scale variations, PDFs, as well as the uncertainties from the non-perturbative corrections, 
which have all been added in quadrature. 
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are missing in the PYTHIA8 model, the cross-section difference c^lo ~~ ^lo * s added to the 
LO QCD+CI prediction in each Mjj and ^dijet bin. With this procedure, we obtain a QCD+CI 
prediction where the QCD terms are corrected to NLO while the CI terms are calculated at 
LO. Non-perturbative corrections due to hadronization and multiple parton interactions are 
applied to the predictions. In Fig. [TJ the predictions are shown for QCD+CI from Ref. |(T2| 
at the CI scale A^ L / RR = 7TeV for the two highest Mjj ranges. The highest Mjj range is the 
most sensitive to the CI signal, though omitting the second-highest Mjj range would decrease 
the expected limit on the CI scale by about 1% for A^ L/RR and 13% for A^ L/RR . Varying the 
boundary between the two highest Mjj ranges by ±0.2 TeV changes the expected limit by less 
then 0.5%. The predictions for the various QCD+CI models at the scale of A = 7 TeV are 
shown in Fig. |2] for the highest Mjj range. At low ^dijet/ the CI predictions with exact NLO 
QCD corrections show smaller enhancement relative to QCD than the corresponding LO CI 
predictions, as described in detail in Ref. [12\. 
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Figure 2: The normalized dijet angular distribution in the highest dijet mass range M« > 
3 TeV compared to various contact interaction models (dashed and dotted histograms). Non- 
perturbative corrections due to hadronization and multiple parton interactions are applied to 
the predictions. The vertical error bars represent the statistical and systematic uncertainties 
added in quadrature. The horizontal bars correspond to the ^dijet bin width. The shaded band 
represents the NLO QCD prediction and includes the systematic uncertainties due to ]i r and ]i / 
scale variations and PDF uncertainties, as well as the uncertainties from the non-perturbative 
corrections added in quadrature. 

The statistical method used to set limits on A follows a modified frequentist approach |3j |33T - 
|35| . The log-likelihood-ratio q = —2 ln( jf 33 ^ 1 ) is used to discriminate between the QCD-only 
hypothesis and the QCD+CI hypothesis. The Lqcd+ci an d ^qcd are written as a product of 
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Poissonian likelihood functions for each bin in ^dijet and for the two highest ranges of M«, 
where the predictions for each M« range are normalized to the number of observed events in 
that range. The p-values, Pqcd+ci(<? > <7obs) an d Pqcd(<7 < <?obs)/ are obtained from ensem- 
bles of pseudo-experiments for the two hypotheses. Systematic uncertainties are represented 
by nuisance parameters which affect the ^dijet distribution. The nuisance parameters are var- 
ied within their Gaussian uncertainties when generating the distributions of q. The QCD+CI 
model is considered to be excluded at the 95% confidence level (C.L.) based on the quantity 
CL S = Pqcd+ci(<? > <?obs)/(l ~~ Pqcd^ < *?obs)) < 0.05. The observed and expected limits 
at 95% C.L. for the CI models considered are listed in Table [2] and displayed in Fig. |3j All 
the observed limits agree within uncertainties with the expected limits, which are evaluated 
at the median of the test statistics distribution of the QCD-only model. The observed limits 
are slightly higher than the expected limits because, for the range M,v > 3.0 TeV, the measured 
dijet angular distribution at low ^dijet is lower than, although statistically compatible with, the 
QCD prediction. The limits for the CI scale A + L/ , RR are also extracted using an alternative pro- 
cedure in which the data are not corrected for detector effects and instead the MC predictions 
are convoluted with the detector resolutions. The limits obtained are found to agree with the 
quoted ones within 1.5%. 

Table 2: Observed and expected lower limits at 95% confidence level for the contact interaction 
scale A for several quark CI models. 



CI model 


Observed limit (TeV) 


Expected limit (TeV) 


NLO A+ L/RR 


7.5 


7.0^ 


-0.4 
-0.6 


NLOA Il/rr 


10.5 


9.7^ 


-1.0 
-1.7 


A LL/RR 


8.4 


7 q+0.5 


A LL/RR 


11.7 


10.9 


+1.7 
-2.4 


LO A yv/AA 


10.4 


9.5^ 


-0.5 
-1.0 


LO A VV/AA 


14.5 


13.7 


+2.9 
-2.6 


LOA fv-A) 


8.0 


7.8^ 


-1.0 
-1.1 



In summary, normalized dijet angular distributions have been measured with the CMS detector 
over a wide range of dijet invariant masses. The distributions are found to be in agreement 
with predictions of QCD and are used to set lower limits on the contact interaction scale for a 
variety of quark compositeness models, including models with NLO QCD corrections. The 95% 
confidence level lower limits for the contact interaction scale A are in the range 7.5 — 14.5 TeV. 
These results represent the most comprehensive set of limits on the contact interaction scale to 
date. 
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Figure 3: Observed (solid lines) and expected (dashed lines) 95% C.L. limits for the contact 
interaction scale A for the different CI models. The dark (light) gray bands indicate the ±lc 
(±2c) uncertainties on the expected limits. 
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